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Molybdenum supported thin films of ScN, MgO and NaF with a Cu adatom have been studied in the frame¬ 
work of density functional theory. We have observed a charge transfer from the metal/film interface to the Cu 
atom and investigated its relation to surface and interface deformations. We find that a weak interaction between 
the metal and the film is a promising prerequisite for adatom charging. The detailed study of Cu/NaF/Mo and 
NaF/Mo indicates that the distortion of the NaF film caused by the Cu adsorption has essentially anharmonic 
character, as it is coupled to a strong charge redistribution in the system. 


I. INTRODUCTION 

Modern experimental techniques made it possible to grow 
epitaxial films with atomic precision creating a whole new 
class of complex materials with unique properties where in¬ 
terfaces often play the decisive role^“^. Ultrathin films of in¬ 
sulating materials grown on a metal support are one exam¬ 
ple of such materials^"^. In particular, due to the proximity 
of the insulator/metal interface to the vacuum/insulator inter¬ 
face the adsorption properties of such complex substrates can 
be completely different from those of the surfaces of the cor¬ 
responding insulating materials^"^’^"^'^. This makes complex 
substrates very interesting for applications in the field of het¬ 
erogeneous catalysis^"^. 

In this regard, one of the intriguing properties of metal sup¬ 
ported thin films is a charging of neutral atoms adsorbed on 
them. This phenomenon was first found when the charge of 
an Au atom adsorbed on NaCl supported by Cu was manipu¬ 
lated with scanning tunneling microscope (STM)^^. The gold 
atom was reversibly switched between the two states: stable 
Au^ and charged Au“. This experimental finding was also 
supported by density functional theory (DFT) calculations^^. 

A spontaneous adatom charging was for the first time pre¬ 
dicted in DFT calculations for an Au atom deposited onto 
MgO/Mo^. The gold atom on this substrate appeared to be 
negatively charged and its adsorption energy turned out to be 
significantly increased as compared to the adsorption of Au 
on MgO^. These theoretical findings were later supported by 
experimental resultsMore recently the charging effect was 
intensively studied theoretically and experimentally for differ¬ 
ent combinations of metal support and thin film materials, and 
adsorbates (see, for example, Ref. 10, 11, 13, 14, 16-35). 

Summarizing the current knowledge about metal supported 
thin films and their adsorption properties we notice the fol¬ 
lowing facts. Thin films can substantially reduce the work 
function of support metals which should lead to eas¬ 

ier charge transfer. However, the charging of adatoms can also 
happen despite the increase of work function^^. For charging 
adatom should have high enough electron affinity^’ 

The origin of adatom charge was found to be the metal and 
oxide interface as shown for Au/MgO/Ag^^, N02/Mg(9/Ag^^, 
NO 2 /BaO/Pt^^, and N02M/2(93/Ag“^. The electron abstrac¬ 
tion from the interface modifies the oxide/metal interaction 


and increases the adhesion between film and metal^^’^^. It was 
also established that the charging of adatom is a long-range 
phenomenon and could occur for relatively thick films 
Honkala et al. 13 showed that the charge of Au on MgO/Mo 
was almost independent of the thickness of MgO ranging from 
1 to 5 monolayers. 

Although the formation of negatively charged species on 
metal supported ultrathin films is quite a widespread phe¬ 
nomenon, the mechanism behind the charging is still under 
debate. First, it was suggested to be a direct electron tun¬ 
neling from the support metal conduction band states to the 
adatom empty states^’^^’^^’^^. However, many authors avoid 
this concept (see, e.g. Refs. 11, 13, 14, 34, and 36), empha¬ 
sizing instead the role of electrostatic interaction and system 
polarization^ We also notice that charge redistribution 

at metal/insulator interfaces was intensively studied in semi¬ 
conductor physics^^""^^ and a number of useful concepts, such 
as band bending and space charge layer formation were devel- 
oped^^“^^ 

It is known that charging is accompanied by surface distor¬ 
tions around the adsorbed adatom^’Moreover, calculations 
show that there is no charging on “frozen” surfacesThe 
importance of rumpling and surface distortion was most con¬ 
sistently studied by Goniakowski et al. (Refs. 34 and 35), who 
showed that the charge transfer from the metal into the thin 
film or adatom is quasilinearly proportional to the rumpling 
of the interface. Moreover, these works reported an important 
finding that the relationship between the charge transfer and 
film deformation were similar for both bare supported films 
and the films with adatoms^^’^^ . Their results suggest that in¬ 
terface rumpling is a response of the system to the sponta¬ 
neous charge redistribution at the interface^"^’^^ . 

Here we report the results of a systematic study of the 
ScN/Mo, MgO/Mo and NaF/Mo systems with a Cu adatom. 
MgO, ScN and NaF were chosen for their simple rock-salt 
structure and small lattice mismatch with Mo. The anions of 
these compounds are neighbours in the Periodic Table of ele¬ 
ments that can yield a trend in properties. Moreover, NaF and 
MgO are ionic compounds, while ScN is more covalent and is 
a semiconductor, which allows us to do a comparative study. 
The questions we focus on are why and when charging takes 
place, how it happens, which atoms and electronic states are 
involved in the charging process. 
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This article is organized as follows. Section II describes 
methodology and computational details. Section III contains 
results and discussions. In Section III A and IIIB we describe 
systems AB/Mo and Cu/AB/Mo (AB = ScN, MgO, NaF), re¬ 
spectively, in particular, factors enhancing charging as well as 
connection between the charge transfer and deformation of the 
film induced by adsorption. In Section III C we investigate the 
origin of the adatom charge and where it accumulates. Next, 
in Section III D we demonstrate how the charge is “pumped” 
from film/metal interface to the top surface layer, and show 
which electronic states are involved. In Section IV we sum¬ 
marize our findings and provide concluding remarks. 


II. METHODOLOGY AND COMPUTATIONAL DETAILS 
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A. Model Description 

The ultrathin films of three binary compounds on the sur¬ 
face of bcc molybdenum (AB/Mo, AB = ScN, MgO, NaF) 
both bare and with an adsorbed Cu atom (Cu/AB/Mo) have 
been studied. Such epitaxially grown thin films usually adopt 
rocksalt structure'^^’'^^, which we consider here. Thus, an 
AB/Mo interface is constructed of the AB (001) surface and 
Mo (001) surface, rotated by 45° with respect to each other. 
The anion atoms (B = N, O, F) are located directly above the 
Mo atoms, while the cation atoms (A = Sc, Mg, Na) are above 
the hollow sites of the Mo lattice. The calculations were car¬ 
ried out for the 2x2 symmetric slabs, having two AB mono- 
layers (ML) on both sides of a 9 ML thick molybdenum slab 
(Fig. 1). The thickness of the Mo slab was chosen based on 
the convergence of the following characteristics i) distances 
between adjacent (001) Mo planes, ii) density of states (DOS) 
of surface Mo layer, in) surface energy. Also, the middle layer 
of the slab should possess bulk-like properties. Our tests show 
that to have a good representation of the mentioned properties 
one can use 9 ML thick Mo slabs. The repeated slabs were 
separated by at least 17 A of vacuum. 

Since the films are ultrathin, while the metal support is sup¬ 
posed to be bulk-like^^ we matched the lattice parameters of 
the films to the one of Mo and aligned them with V2 x <20 (Mo) 
(^^2 factor appears due to the Mo(OOl) surface rotation by 
45°). The lattice mismatch parameters = af/as - I (af 
and as are lattice constants of the film and support) for ScN, 
MgO and NaF on Mo are 1.2%, -5.3% and 4.1%, respec¬ 
tively (as calculated from the experimental lattice constants 
Refs. 44-47). Plus sign here means that the film is com¬ 
pressed with respect to the bulk lattice parameter, while minus 
means the film is expanded. The lattice mismatch is modest 
and the most strained film, MgO/Mo, has successfully been 
synthesized (see, e.g. Ref. 6 and 43). 

For all the surfaces we considered four sites for Cu adsorp¬ 
tion: 1) on top of cation, 2) on top of anion, 3) hollow site 
above the center of the square comprised of two cations and 
two anions, and 4) bridging position above the center of the 
cation and anion bond (Fig. 1). The Cu atom at these sites 
was allowed to relax along the Z axis. 


FIG. 1: (Color online) Left panel: View of the upper half of a 
2x2 symmetric Cu/MgO(2ML)/Mo(9ML)/MgO(2ML)/Cu 
slab. Monolayers are enumerated with index j. The dashed 
line represents the MgO/Mo interface and the middle 
symmetry layer of Mo has index 7 = -5. The vacuum gap 
between repeated slabs is not shown. Right panel: Top view 
of the MgO surface showing four adsorption sites: 1) Mg 
on-top, 2) O on-top, 3) hollow, 4) bridging. 


B. Computational Details 

Density functional theory (DFT) calculations were per¬ 
formed using the projector augmented wave method"^^ to¬ 
gether with the general gradient approximation (GGA) 
in Perdew-Budke-Emzemhof parametrization'^^’^^ as imple¬ 
mented in Vienna Ab initio Simulation Package^ The 
cutoff energy of 800 eV was used in the calculations. The 
following states were included into the valence band: Cu 
3Ji04^i, Sc 3JM^^ N 2s^2p\ Mg 3^^ O 2s^2p\ Na 2 / 3 ^ 1 , 
F 2s^2p^, and Mo . The Brillouin zone was sampled 

over a 8 X 8 X 1 Pack-Monkhorst k-points mesh^^. The calcu¬ 
lations were spin polarized. During geometry optimization the 
relaxation of all atoms except for the Mo atoms in the middle 
symmetry layer (7 = -5, see Fig. 1) was allowed. A geometry 
was considered to be optimized if forces acting on unfrozen 
atoms were less than 5 meV/A. 

We notice that even though standard GGA functionals un¬ 
derestimate band gaps, they adequately describe charging 
phenomena in metal supported thin films that has been demon¬ 
strated by numerous previous More¬ 

over, it has recently been shown that both standard GGA and 
hybrid functional produce very similar Bader charges of an 
Au adatom on MgO and CaO substrates doped with Mo or 

To understand the mechanism of adatom charging on the 
AB/Mo surfaces we analysed the atomic charges calculated 
using the Bader approach^^’^^. We also found it elucidative 
to analyze the sum of atomic charges in separate layers (an 
approach also used in Ref. 60). In order to choose an optimal 
grid we tested the convergence of Bader charges with respect 
to the grid density. Our tests showed that the Bader charge of 
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a monolayer was much more sensitive to the grid step in the 
Z direction than in the other directions that allowed us to save 
computational resources by choosing a non-uniform grid with 
smaller steps along the Z axis. The step values: Az < 0.02 A 
and Ax = Ay < 0.05 A resulted in the Bader charges with the 
accuracy of about 0.01 - 0.05^ per monolayer and about 0.01^ 
per atom. 


III. RESULTS AND DISCUSSIONS 
A. ScN/Mo, MgO/Mo, NaF/Mo 

First, we calculated the bulk and surface properties of the 
involved materials. The calculated lattice constant of bcc Mo, 
gq = 3.151 A, is in good agreement with the experimen¬ 
tal value of 3.147A'^'^, and the calculated surface energy of 
Mo(OOl), 3260mJ/m^ (0.203 eV/A^), is close to experimen¬ 
tal 2930 mJ/m^ (0.183 eV/A^)^^ The calculated work func¬ 
tion of Mo(OOl), O = 3.76 eV, is underestimated compared 
to the experimental value of 4.53 eV^^^. The underestimation 
of the work function by GGA for similar systems is well- 
known^^. The calculated lattice constants of bulk ScN, MgO 
and NaF am 4.500 A (4.501 A^^), 4.240 A (4.214 A^^), 4.705 
A (4.634 respectively (experimental values are given in 
parentheses). Notice that in the calculations of the complex 
substrates reported here the Mo lattice parameter was used. 
The calculated band gap values for ScN, MgO and NaF are 
0.1 eV (1.3 eV^^3), 4.6 eV (7.7 eV^^) 5 ^ 5 

respectively (experimental values are given in parentheses). 
As expected the band gap values are underestimated. 

The electronic structure, charge density redistribution pat¬ 
tern and other calculated characteristics of the metal supported 
thin films are shown in Figs. 2, 3 and Table I. The densi¬ 
ties of states demonstrate that the states of ScN and Mo are 
strongly mixed, the states of MgO and Mo are moderately 
mixed, while for those of NaF and Mo virtually there is no 
mixing (Fig. 2). The calculated adhesion energy follows this 
trend decreasing from ScN to NaF. The interface distance, on 
the contrary, increases from ScN to NaF (Table I), indicating 
a somewhat stronger bonding between ScN and Mo compared 
to the other considered cases. The charge density difference 
maps (Fig. 3) show a noticeable charge redistribution at the 
ScN/Mo interface involving even the second layer of ScN. In 
the case of MgO/Mo the charge redistribution at the interface 
is less pronounced with mostly O atoms being affected. In the 
case of NaF/Mo the redistribution is quite weak. Rumpling at 
the interface is more pronounced for MgO/Mo and NaF/Mo, 
than for ScN/Mo. The reduction of the work function is least 
pronounced for ScN/Mo and most pronounced for NaF/Mo. 

Thus, our results demonstrate that the strength of Mo and 
film bonding decreases in the row: ScN, MgO, and NaF 
(Fig. 2, Fig. 3, and Table I). Among the considered com¬ 
pounds ScN is the least ionic and NaF is the most ionic, 
which is also shown by the analysis of the Bader charges of 
the bulk compounds: 2(N) = -1.60^, 2(0) = -1.65^, and 
2(F) = -0.87^. Thus we observe the following trend: for 
more ionic films we find weaker adhesion to the metal, smaller 


TABLE I: Adhesion energy L’adh, work function O, interface 
distance d = ^Zo(A) -r ^zo(B) - z_i(Mo), rumpling in the 
layer ry = Zj(B) - Zj(A). 


ScN/MO MgO/Mo NaF/Mo 


£adh (eV/A2) 

0.204 

0.130 

0.039 

OCeV) 

3.14 

1.85 

1.28 

d(A) 

2.11 

2.15 

2.41 

n (y 

0.02 

0.04 

0.02 

'•0 (A) 

0.02 

0.08 

0.04 


overlap of the densities of states of the metal and the film, and 
less pronounced charge redistribution at the interface. As we 
show further, the strength of the interface bonding is also re¬ 
lated to the charging of an adatom on the AB/Mo surface. 


B. Cu Atom Adsorption on Mo Supported Thin Films 

Spontaneous adatom charging is studied here by consider¬ 
ing the adsorption of Cu atom on the AB/Mo surfaces. Be¬ 
fore going into details, we want to outline that the capability 
to negatively charge the Cu adatom increases in the sequence 
ScN/Mo, MgO/Mo, NaF/Mo — the same sequence in which 
the bonding strength, adhesion energy, intensity of charge re¬ 
distribution and the work function of the AB/Mo substrate de¬ 
crease. 

In order to characterize the adsorption of Cu we investi¬ 
gated four adsorption sites: on top of cation, on top of an¬ 
ion, hollow and bridging, and calculated several quantities, 
such as, adsorption energy F’ads, Bader charge of Cu adatom 
2(Cu), the vertical distance from Cu to the top layer (av¬ 
erage position of atoms) z, and rumpling of atomic layers 
C “ - (Z;(A)), calculated from average z coordinates 

of A and B in the layer with an index j (see Fig. 1). The 
obtained results are summarized in Table II. 

As one can see charging of Cu is most pronounced on 
NaF/Mo. For all the considered adsorption positions, the Cu 
4s states fall below the Fermi level and no spin-polarization 
is observed (see the example of F on-top site. Fig. 4c). The 
Bader charge of Cu is almost the same for all the positions, 
2(Cu) -0.16e, which is in agreement with charges reported 

for similar systems, e.g. Au on MgO/Mo^^. Interestingly, de¬ 
spite fluorine being the most electronegative element of the 
Periodic Table Cu charging occurs even for Cu adsorbed on- 
top of F, which was found to be the most favorable adsorption 
site. 

Charging of Cu on NaF/Mo is accompanied by substantial 
atomic displacements in the top NaF layer that results in rather 
large adsorption energies (from 1.651 to 1.877 eV). In partic¬ 
ular, the F atom, situated below adsorbed Cu, moves down¬ 
wards by 0.27 A, while the four nearest Na atoms move up¬ 
wards by 0.57 A. A similar distortion pattern was reported 
for Au/MgO/Mo^. Together with distortions around Cu we 
also observe a prominent average rumpling -0.34 A) of 
the NaF surface layer (j = 1, see Fig. 1) for all the adsorp- 
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FIG. 2: (Color online) The partial density of states of (a) ScN/Mo, (b) MgO/Mo, and (c) NaF/Mo. The Fermi level is set to zero. 
From left (ScN/Mo) to right (NaF/Mo) the overlap of metal and film states decreases drastically. For Sc total DOS (T) is shown. 
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(a) ScN/Mo (b) MgO/Mo (c) NaF/Mo 


FIG. 3: (Color online). Charge density difference maps of Ap = p(AB/Mo) -p(AB) -p(Mo) shown for (100) plane: (a) 
ScN/Mo, (b) MgO/Mo, (c) NaF/Mo. Here, red represents an increase and blue a decrease of charge density. The intensity of 
charge redistribution decreases from left (ScN/Mo) to right (NaF/Mo). 


tion sites (see Table II). The significant changes occur also at 
the interface, where the inversion of rumpling is observed and 
Zo(F) - z-i (Mo) becomes smaller than zo(Na) - z-i (Mo). Thus 
upon Cu adsorption on-top of F the Mo-F distance on average 
shortens by ~ 0.1 A, as compared to that in NaF/Mo. For all 
the adsorption sites on NaF/Mo the same inversion of rum¬ 
pling occurs (see Table II). Similar structural changes were 
also reported for the NO 2 adsorption on A\ 20 ^IAg^^. 

The characteristics of Cu adsorption on MgO/Mo are sim¬ 
ilar to those of Cu/NaF/Mo. Cu charging occurs for all the 
four adsorption sites. The Bader charges of Cu, 2(Cu), are 
smaller than those for Cu on NaF/Mo, but still significant 
ranging from -0.59^ to -Q.61e. The adsorption energies are 


rather large (from 1.326 eV to 1.466 eV), the favored adsorp¬ 
tion site is the hollow site. The Cu adsorption on MgO/Mo 
is accompanied by notable surface relaxation. The two Mg 
atoms, nearest to Cu adsorbed at the hollow site, shift up¬ 
wards by 0.19 A, while two nearest O atoms move downwards 
by 0.11 A. The average rumpling of the surface monolayer 
(7 = 1) is about -0.07 A. At the same time, rumpling of 
the interface MgO layer (7 = 0) is negligible as compared to 
that in the Cu/NaF/Mo system. The average vertical distances 
from O (7 = 0) to the Mo (7 = -1) layer and from Mg (7 = 0) 
to the Mo (7 = -1) layer are the same, 2.15 A. Due to the Cu 
adsorption the Mo-0 distance shortens by ~ 0.04 A. 

According to Bader analysis the Cu atom on MgO/Mo has 





























































5 



FIG. 4: (Color online). The density of states of (a) N on-top site Cu/ScN/Mo, (b) hollow site Cu/MgO/Mo, and (c) F on-top site 
Cu/NaF/Mo. The Fermi level is set to zero. 


TABLE II: Cu adsorption on AB/Mo (AB = ScN, MgO, 

NaF) considering four different adsorption positions. 
Adsorption energy £’ads, Bader charge of Cu adatom g(Cu), 
vertical distance z from Cu to the top layer (average position 
of atoms). The average rumpling is calculated from the 
average z coordinates of A and B in the layer. The adsorption 
on frozen surfaces (only Cu atom position was optimized) is 
marked with stars (*). Bridging position on the ScN/Mo 
substrate was found to be unstable and was not converged. 


£ads (eV) e(Cu) (e) Z (A) ro (A) n (A) 

ScN/Mo 

Sc on-top 

0.752 

-0.51 

2.84 

-0.02 -0.06 

N on-top 

1.609 

-0.01 

1.96 

0.04 0.03 

N on-top* 

1.496 

-0.07 

1.93 

0.02 0.02 

hollow 

0.909 

-0.39 

2.36 

-0.01 -0.03 

MgO/Mo 

Mg on-top 

1.386 

-0.64 

2.77 

0.00 -0.07 

0 on-top 

1.326 

-0.59 

2.26 

0.01 -0.07 

hollow 

1.466 

-0.67 

2.45 

0.00 -0.07 

hollow* 

0.848 

-0.32 

2.30 

0.08 0.04 

bridging 

1.401 

-0.65 

2.50 

0.00 -0.07 

NaF-Mo 

Na on-top 

1.651 

-0.76 

2.78 

-0.22 -0.34 

F on-top 

1.877 

-0.76 

2.69 

-0.22 -0.34 

F on-top* 

0.509 

-0.22 

2.47 

0.04 0.02 

hollow 

1.740 

-0.76 

2.92 

-0.21 -0.34 

bridging 

1.764 

-0.76 

2.80 

-0.21 -0.34 


a smaller charge than that on NaF/Mo, however, in both cases 
the Cu 4s states are filled and no spin-polarization is observed 
(Fig. 4). Although the Bader analysis cannot be expected to 
unambiguously assign charges to atoms in a solid somewhat 
smaller charge of Cu on MgO/Mo corresponds to a smaller 
adsorption energy and smaller displacements of surface atoms 
in comparison with those for NaF/Mo. 

On ScN/Mo, Cu prefers to adsorb on top of N that yields no 
charging. The Cu 4^' states are half occupied and mixed with 
the N states (Fig. 4a). Adsorption at this site results in Cu-N 


covalent bonding with the Cu-N distance of 1.86 A and the 
adsorption energy of 1.609 eV. Surface relaxation in this case 
is very weak in comparison with all the cases where charging 
is observed. 

As a matter of fact, charging of Cu occurs also on ScN/Mo 
for two adsorption sites: Sc on-top and hollow. However, even 
in these cases it is less pronounced than in the Cu/MgO/Mo 
and Cu/NaF/Mo systems. For Cu adsorbed at the Sc on- 
top site the Cu 4^’ states are fully occupied (not shown) and 
at the hollow site Cu 4^' states are preponderantly occupied 
(not shown). The calculated Bader charges are still notable, 
2(Cu) = -0.51^ and -0.39^ for the Sc on-top and hollow 
sites, respectively. The interface and surface rumpling is also 
negative as in the case of Cu/NaF/Mo but somewhat less pro¬ 
nounced (see Table II). 

We notice that as Cu gets charged at all adsorption sites on 
both MgO/Mo and NaF/Mo the resulting adsorption energies 
are quite similar that provides a fiat energy profile for Cu on 
these surfaces, especially on MgO/Mo, and, therefore, simpli¬ 
fies surface diffusion. This is not the case for ScN/Mo where 
the adsorption energies are very different, e.g. 0.752 eV and 
1.609 eV for the Sc on-top and N on-top positions, respec¬ 
tively. 

It is well-known that without metal support adatom charg¬ 
ing is negligible (see, e.g. Ref. 9 and 66). We have examined 
Cu adsorption on top of anions in the Cu/AB systems (2ML 
AB with the same lattice parameter as in the other calcula¬ 
tions). For the N on-top site of Cu/ScN we got £’ads = 1.523 
eV and 2(Cu) = -0.13^. For the O on-top site of Cu/MgO we 
got L’ads = 1.110 eV and g(Cu) = -0.16^. For the F on-top 
site of Cu/NaF we got ^'ads = 0.380 eV and 2(Cu) = -0.12^. 
Thus, in all these cases charging is negligible. 

The comparison of Cu on ScN/Mo, MgO/Mo and NaF/Mo 
allows us to make some general conclusions regarding the 
capability of metal supported thin films to charge adatoms. 
As follows from our results this capability is enhanced from 
ScN/Mo to MgO/Mo and further to NaF/Mo. ScN is a co¬ 
valent compound also forming covalent-like bonds with Mo, 
while MgO and NaF are ionic (see Section III A). There¬ 
fore, to construct a system where a substantial charge trans- 
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FIG. 5: (Color online) Charge of Cu adatom in Cu/AB/Mo 
versus energy of surface distortion, 

^dist = ^([Cu]/AB/Mo) - ^(AB/Mo). 


fer is possible one could use ionic films whose states are 
weakly mixed with those of the metal support. Strong adatom- 
substrate bonds can also prohibit the charge transfer whereas 
a weaker interaction between the adatom and a film can pro¬ 
mote it. Notice that Cu and MgO (NaF) states almost do not 
mix. The same trend we find for Cu adsorbed at the Sc on- 
top and hollow sites on ScN/Mo where a charge transfer takes 
place. Thus, a minimum mixing between adatom and thin film 
states as well as between those of the insulating film and the 
metal support appear to be crucial for charging to occur. 

We also investigated the adsorption of Cu on frozen 
surfaces for the most stable positions: N on-top site for 
Cu/ScN/Mo, hollow site for Cu/MgO/Mo, and F on-top site 
for Cu/NaF/Mo (see Table II). For the F on-top (Cu/NaF/Mo) 
and hollow sites (Cu/MgO/Mo) the adsorption energies and 
Bader charges of Cu are reduced significantly. According to 
the density of states calculated for these cases the Cu spin- 
down 4s states are split from the spin-up states and are not 
fully occupied. These findings agree with the results reported 
for the Au adsorption on frozen MgC/Ag^"^’^^. We also no¬ 
tice that in the case of Cu adsorption at the N on-top site of 
ScN/Mo the adsorption energy, Bader charge and the Cu-N 
distance calculated for the frozen and fully relaxed surfaces 
are very similar (see Table II). This is due to the absence of 
charging in both cases, when the Cu-N bond is a dominating 
contribution to the adsorption energy. Also we observe, that if 
there is no charging or it is significantly reduced (as on frozen 
surfaces) the vertical distance from the Cu to the surface is 
shorter, while charged Cu shifts further away from the surface 
(see Table II). 

Thus our results confirm the importance of surface relax¬ 
ation for adatom charging. We also find it useful to study 
the characteristics of the AB/Mo systems having exactly the 
same surface distortion as the Cu/AB/Mo systems (i.e. when 
all A, B and Mo atoms are kept in the same positions as in 
Cu/AB/Mo, but Cu is removed). We denote such systems by 


[Cu]/AB/Mo. The energy needed to distort the [Cu]/AB/Mo 
surface is then £’dist = £’([Cu]/AB/Mo) - £’(AB/Mo). £’dist 
could be used as a measure of surface distortion. In Fig. 5 we 
show the charge of Cu in Cu/AB/Mo versus distortion energy 
of [Cu]/AB/Mo, E’dist- There is a clear correlation between the 
energy of distortion and charge transferred to Cu, in fact, £’dist 
and charge are linearly proportional to each other. 


C. Distortion Induced Charge Localization on AB/Mo surface 

The charge accumulated by the Cu adatom comes from 
AB/Mo, but where is the origin of this charge? The anal¬ 
ysis of the Bader charges of the substrate monolayers, cal¬ 
culated as described in Section IIB, makes it possible to an¬ 
swer this question. Here we study the charge distribution in 
Cu/AB/Mo, [Cu]/AB/Mo and AB/Mo. In Fig. 6 we show the 
layer charges for the three most stable configurations: N on- 
top site for Cu/ScN/Mo, hollow site for Cu/MgO/Mo and F 
on-top site for Cu/NaF/Mo. 

Let us, first, compare Cu/AB/Mo and AB/Mo. Fig. 6 (top 
panels) shows that for all Cu/AB/Mo the charge redistribu¬ 
tion mostly occurs at the interface, while deeper Mo layers 
remain practically unaffected. For Cu adsorbed at the N on- 
top site, neither significant changes at the interface can be 
seen. For Cu adsorbed at the hollow site on MgO/Mo, where 
a significant charge transfer to Cu takes place, the charge of 
the adatom originates from the metal/oxide interface, mostly 
from the first Mo layer (Fig. 6b top panel). These findings 
for Cu/MgO/Mo and Cu/NaF/Mo agree well with previous re¬ 
ports on Au/MgO/Mo and N02/MgO/Mo^'^. 

Let us now consider distorted [Cu]/AB/Mo, which has 
the geometry of relaxed Cu/AB/Mo but without Cu. In 
terms of charge redistribution not much difference between 
Cu/ScN/Mo and [Cu]/ScN/Mo is observed (see Fig. 6a). For 
[Cu]/MgO/Mo and [Cu]/NaF/Mo, however, one can see a 
strong polarization of AB and Mo as a result of the distor¬ 
tion. Similar to Cu/MgO/Mo and Cu/NaF/Mo, the charge is 
accumulated in the top surface monolayer of MgO and NaF 
(see Fig. 6). The amount of charge accumulated in the surface 
layer is -0.39^ for MgO and -0.77e for NaF. The origin of the 
charge is again the MgO/Mo (NaF/Mo) interface (see Fig. 6). 

The Bader analysis of atomic charges shows that the excess 
charge, which for a particular atom we define as 2[Cu]/ab/Mo “ 
2 ab/Mo5 is spread over several atoms in the top layer (Fig. 6, 
bottom panels). In [Cu]/MgO/Mo most of the excess charge is 
accumulated in the surface layer around the hollow site and in 
the next coordination sphere. In [Cu]/NaF/Mo -0.32^ excess 
charge is localized at the F atom (below the Cu in Cu/NaF/Mo) 
and -0.26^ is localized at the four surrounding Na atoms. 
Thus, for both [Cu]/MgO/Mo and [Cu]/NaF/Mo the excess 
charge appears to be spread over the distorted top layer and 
its larger part is located at the adsorption site. The amount 
of this excess charge accumulated in the distorted AB top 
layer is almost equal to the charge of Cu in the corresponding 
Cu/AB/Mo system. Thus, in both cases with and without Cu, 
the distortion leads to similar charge redistribution. For more 
details on charge redistribution see Fig. 8 in Appendix A. 
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FIG. 6 : (Color online) Top panels: Comparison of charges of different layers of AB/Mo, Cu/AB/Mo and [Cu]/AB/Mo (has 
frozen geometry of relaxed Cu/AB/Mo, but Cu is absent). Bottom panels: Accumulation of excess charge on specific atoms (or 
groups of atoms) in the top layer of AB, i.e. 2[Cu]/ab/Mo “ 2ab/Mo- (a) N on-top site Cu/ScN/Mo, (b) hollow site Cu/MgO/Mo, 
and (c) F on-top site Cu/NaF/Mo. 


D. Charge Transfer From Metal Support to Adatom Through 
Thin Film: Cu/NaF/Mo case 

In this section we discuss the role of the thin film distortion 
in the charge transfer to the Cu adatom using the case of Cu 
adsorbed at the F on-top site on NaF/Mo (Cu/NaF/Mo), where 
charging is most pronounced. 

To see how charge redistribution depends on the degree of 
deformation we approximated the deformation path of the sur¬ 
face with a linear geometry transformation. Let us denote 
the set of coordinates of all the Na and F atoms in NaF/Mo 
with Ro and those in Cu/NaF/Mo with Ri, and introduce 
Ai^ = Ri - Rq. Next, let us consider the linear path of 
geometry transformation from R^ to Ri, which is given by 
R = Ro + ^AR, where 0 < ^ < 1. For every value of ^ the 
positions of Na and F atoms are defined and frozen, while Cu 
is allowed to relax. 

Furthermore, we have considered the same linear transfor¬ 
mation of NaF/Mo in the absence of Cu (the [Cu]/NaF/Mo 
case), as well as that of the NaF film in the absence of both 
Cu and Mo (the [Cu]/NaF/[Mo] case). Results are shown in 
Fig. 7. 

It is clear that in the case of Cu adsorption on NaF/Mo 
(Cu/NaF/Mo) the relaxation occurs spontaneously without 


any barrier. The charge of adatom increases almost linearly 
along the deformation path. Bader analysis shows that for 
both Cu/NaF/Mo and [Cu]/NaF/Mo the charges of NaF (j = 0, 
see Fig. 1) and Mo (j = 1) interface layers monotonically de¬ 
crease along the path (shown for [Cu]/NaF/Mo, see Fig. 7b), 
thus charge is depleted from the NaF/Mo interface. 

For Cu/NaF/Mo the gradually applied distortion yields 
gradual charging of Cu, while in the case of [Cu]/NaF/Mo it 
provides monotonic charge accumulation in the NaF top layer 
(j = 1 ) (see Fig. 7b). Moreover, we notice that along the path 
the values of the charge of Cu in Cu/NaF/Mo and the charge of 
NaF (7 = 1) top layer in [Cu]/NaF/Mo are similar and become 
almost equal for ^ = 1 . 

Furthermore, our results indicate that the Nsl 3s and ¥ 2p 
states participate in the charge transfer from the NaF/Mo in¬ 
terface to the Cu adatom (see Fig. 9 in Appendix B). 

Our calculations show that charging requires distortion, 
i.e. cannot happen immediately (see Fig. 7). There are in¬ 
dications that metal supported ultrathin films have different 
phonon structure compared to thick films and should be more 
fiexible^^. Indeed, energy changes shown in Fig. 7a demon¬ 
strate that it is much easier to deform [Cu]/NaF/Mo than 
[Cu]/NaF/[Mo]. The same distortion from^ = 0 to ^ = 1 costs 
0.70 eV for [Cu]/NaF/Mo and 2.58 eV for [Cu]/NaF/[Mo]. 
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FIG. 7: (Color online) Linear geometry transformation path from ^ = 0 to ^ = 1 : (a) energy change along the path for 
Cu/NaF/Mo, [Cu]/NaF/Mo and [Cu]/NaF/[Mo], (b) charge of Cu adatom in Cu/NaF/Mo, and charges of NaF (j = 1) top 
surface layer, NaF (j = 0) interface layer and Mo (j = -1) interface layer in [Cu]/NaF/Mo (in the absence of Cu). 


Moreover, the analysis of the curves show rather harmonic 
behavior for [Cu]/NaF/[Mo], while for [Cu]/NaF/Mo and 
Cu/NaF/Mo it is strongly anharmonic (see Fig. 7a). Notice 
that with (or without) Cu adatom the energy change is almost 
linearly proportional to the accumulated charge. Hence, we 
attribute the NaF/Mo anharmonicity to the coupling of de¬ 
formation with the charge transfer from the Mo/NaF inter¬ 
face to the NaF top layer (in [Cu]/NaF/Mo) or Cu adatom (in 
Cu/NaF/Mo). 


IV. CONCLUSIONS 

We have performed a systematic study of the AB/Mo and 
Cu/AB/Mo systems (AB = ScN, MgO, NaF). An enhance¬ 
ment of Cu adatom charging from ScN to MgO, and further 
to NaF is observed. The results suggest that charging is more 
pronounced when mixing between the states of the insulating 
film and metal substrate as well as between those of an adatom 
and thin film is small. 

The film/metal interface is the origin of adatom charge. 
The charge transfer is accompanied with a strong surface 
relaxation around the adatom and structural changes at the 
film/metal interface. Our results on Cu/AB/Mo show clear 
correlation between the amount of transferred charge and the 
degree of system distortion. Moreover, the results of the con¬ 
strained deformation calculations indicate that distortion is 
decisive for the charge transfer and it will “pump” the charge 
into the top surface layer irrespective of an adatom presence 
or absence on the surface. 

We have found that the deformation of NaF/Mo film is 
essentially anharmonic, which we attribute to the coupling 
with the charge “pumping” from metal/film interface to the 
film/vacuum interface or adatom. Also, it is much easier to 
distort NaF/Mo than NaF, therefore we expect that soften¬ 


ing of some phonon modes of the metal supported films takes 
place. 
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Appendix A: Spatial Charge Distribution 

The spatial distribution of charge accumulated in 
[Cu]/NaF/Mo was examined further. In Fig. 8 we show 
charge density for [Cu]/NaF/Mo and [Cu]/NaF/[Mo] inte¬ 
grated along the Z axis. Let us recall, that [Cu]/NaF/Mo (also 
[Cu]/NaF/[Mo]) have the same geometry as Cu/NaF/Mo, 
but Cu (respectively Cu and Mo) are absent. One can 
see that in the [Cu]/NaF/Mo case charge is more spread 
in the Z direction compared to the case when both Cu 
adatom and Mo support are removed ([Cu]/NaF/[Mo]). In 
particular, for [Cu]/NaF/Mo approximately 0.2e are found 
above the position of Cu in Cu/AB/Mo (note, there is no Cu 
in [Cu]/NaF/Mo). The same conclusion follows from the 
integrated magnetization curve, which agrees well with the 
integrated charge curve (Fig. 8). 
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FIG. 8: The comparison of [Cu]/NaF/[Mo] with 
[Cu]/NaF/Mo showing charge spillover for the latter. Charge 
or magnetization integrated between z-Z' plane and the 
middle of the vacuum gap. Z' = 0 corresponds to the z 
coordinate of [Cu]/NaF/Mo surface termination. 


Appendix B: DOS of Na and F during deformation 

Here we report the details behind the distortion of 
[Cu]/NaF/Mo from ^ = 0 to ^ = 1 (see the main text for ex¬ 
planations). In Fig. 9 we show the evolution of the density of 
states near the Fermi level for [Cu]/NaF/Mo, namely, for the F 
atom underneath Cu in Cu/NaF/Mo and the nearest Na atoms 
in the top NaF layer. According to the Bader analysis these 
atoms accumulate the charge, while Mo atoms in the interface 
layer lose it (not shown) that agrees well with the evolution 
of the density of states. The number of the F states below 
the Fermi level (mainly the states) and Na states (mainly 
the 3 5' states) gradually increases along the deformation path 
(see Fig. 9), while the number of the Mo states decreases (not 
shown). As a result of the charge transfer the system becomes 
spin polarized, near the Fermi level the spin-down states are 
unoccupied (see Fig. 9). Therefore, our findings indicate that 
the F 2/7 and Na states participate in the charge transfer to 
Cu adatom. 
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FIG. 9: [Cu]/NaF/Mo: partial density of states near the Fermi level (set to zero): (a), (b), and (c) Na atom (which is near Cu 
adatom in Cu/NaF/Mo), (d), (e), and (f) F atom (which is below Cu adatom in Cu/NaF/Mo). In (a) and (d) ^ = 0.0, in (b) and 
(e) ^ = 0.5, in (c) and (f) ^ = 1.0. 
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